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Abstract 

Oscillating  solutions  to  the  effective  equations  of  Loop  Quantum  Cosmology  have  been  suggested 
for  the  role  of  an  ‘eternal  seed’,  providing  a  possible  starting  point  for  the  emergent  universe  sce¬ 
nario.  We  investigate  the  stability  of  a  particular  model,  sourced  by  a  homogeneous  massless  scalar 
field  and  a  negative  cosmological  constant,  with  respect  to  small  perturbations  and  to  quantum 
collapse.  We  find  that  the  model  has  perturbatively  stable  and  unstable  solutions,  with  both  types 
of  solutions  occupying  significant  regions  of  the  parameter  space.  All  solutions  are  unstable  with 
respect  to  collapse  by  quantum  tunneling  to  zero  size.  We  discuss  the  possibility  that  the  state  re¬ 
sulting  from  the  collapse  is  non-singular,  so  it  may  tunnel  back  to  the  oscillating  regime.  We  argue 
that  the  universe  is  then  likely  to  evolve  to  states  of  very  large  size  with  large  particle  occupation 
numbers. 
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I.  INTRODUCTION 


Cosmic  inflation  may  be  eternal  to  the  future,  but  it  cannot  be  eternal  to  the  past.  More 
specifically,  if  the  average  expansion  rate  Hav  along  some  geodesic  is  positive,  that  geodesic 
must  be  past-incomplete  [T] .  This  suggests  that  the  inflationary  phase  of  our  universe  had 
some  sort  of  a  beginning,  perhaps  a  quantum  nucleation  from  ‘nothing’  [2]. 

A  possible  alternative  is  a  universe  which  existed  as  an  eternal  “seed”  before  sprouting  into 
our  inflating  universe.  In  Refs,  p]  and  [1] ,  Ellis  and  collaborators  proposed  a  mechanism  for 
this  type  of  “emergent”  universe  scenario.  Their  model  relies  on  a  scalar  held  in  a  potential 
which  is  asymptotically  hat  as  0  — >  —  oo,  then  reaches  a  hill  as  0  — >  0,  signaling  the  beginning 
of  inflation.  So  long  as  the  condition  Hav  =  0  is  satished  prior  to  the  beginning  of  inhation, 
the  conclusions  of  [Ij  do  not  apply. 

However,  for  the  emergent  universe  to  be  a  viable  scenario,  it  must  be  possible  to  construct 
a  model  which  is  stable  for  eternity.  One  candidate  model  was  discussed  by  Graham  et  al 
[5],  who  found  that  a  spherical  (K  =  +1)  universe  containing  matter  with  equation  of  state 
w  =  —2/3  and  a  negative  cosmological  constant  A  <  0  oscillates  as  a  simple  harmonic 
oscillator.  The  negative  A  causes  an  expanding  universe  to  turn  around  and  recollapse, 
while  the  matter  component  with  w  =  —2/3  causes  a  collapsing  universe  to  bounce.  The 
result  is  oscillation.  A  perfect  fluid  with  w  =  —2/3  would  be  unstable  with  respect  to  high- 
frequency  perturbations,  but  Graham  et  al  argue  that  a  frustrated  network  of  domain  walls 
would  have  an  effective  equation  of  state  w  =  —2/3  and  be  perturbatively  stable.  However, 
it  has  been  shown  in  [6]  that  this  ‘simple  harmonic  universe’  can  tunnel  from  the  bounce 
point  to  nothing,  signaling  a  quantum  instability  and  indicating  that  a  universe  of  this  type 
cannot  last  forever. 

A  different  class  of  models  for  an  emergent  universe  are  motivated  by  Loop  Quantum 
Cosmology  (LQC),  which  provides  a  nonsingular  bounce  scenario  due  to  a  modification  in 
the  effective  Friedmann  equation.  Mulryne  et  al  [T]  have  shown  that  such  a  model  has 
eternally  oscillating  solutions  with  positive  spatial  curvature  ( K  =  +1),  and  that  inflation 
may  arise  with  an  appropriate  potential  for  the  scalar  held.  An  LQC  model  with  hat  spatial 
geometry  ( K  =  0)  and  the  energy  density  given  by  a  massless  scalar  held  0  and  a  negative 
cosmological  constant  was  recently  discussed  by  Miclczarek  et  al  [8j,  who  found  a  simple 
oscillating  solution  in  this  case.  It  should  be  noted  that  Loop  Quantum  Gravity,  on  which 
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LQC  is  based,  is  still  an  incomplete  theory;  in  this  sense  the  foundations  of  LQC  are  not  very 
reliable.  However,  LQC  has  now  been  studied  in  great  detail  and  is  an  interesting  (albeit 
restricted)  theory  in  its  own  right  (see  Ref.  [2]  for  an  up  to  date  review).  It  resolves  the 
singularities  of  FRW  models  and  provides  a  useful  framework  for  investigating  Planck-scale 
physics. 

In  the  present  paper,  we  shall  investigate  the  stability  of  the  oscillating  universe  of  Ref.  [8] 
with  respect  to  classical  scalar  held  perturbations  and  to  quantum  tunneling.  In  the  next 
Section  we  introduce  the  model  and  its  oscillating  solution.  The  stability  of  this  solution 
with  respect  to  small  perturbations  is  analyzed  in  Section  III.  We  find  that  unstable  modes 
are  present  for  a  significant  part  of  the  parameter  space  of  the  model,  but  there  are  also 
significant  ranges  of  the  parameters  where  the  model  is  perturbatively  stable.  This  situation 
is  similar  to  that  found  for  the  simple  harmonic  universe  in  Ref.  [5] ,  but  there  are  also  some 
significant  differences,  as  we  shall  see. 

The  stability  of  the  oscillating  LQC  model  with  respect  to  quantum  tunneling  will  be 
investigated  in  Section  IV.  We  earlier  found  in  Ref.  [10J  that  in  order  to  stabilize  an  oscillating 
universe  against  quantum  collapse  in  standard  Einstein  gravity,  one  would  need  a  matter 
component  with  a  negative  energy  density  that  grows  faster  than  a~6  at  a  — >  0,  where  a  is  the 
scale  factor.  As  we  shall  discuss  in  the  following  Section,  the  effective  Friedmann  equation 
of  LQC  replaces  the  matter  density  p  with  p  =  p(  1  —  p/ pmax ),  where  pmax  =  const.  The 
energy  density  of  a  massless  scalar  held  is  p  oc  a-6,  and  thus  p(a  — *  0)  oc  —  cT12,  suggesting 
the  possibility  that  the  model  could  be  non-perturbatively  stable.  We  find,  however,  that 
the  results  of  Ref.  [TO]  cannot  be  applied  here,  because  the  gravitational  action  in  LQC  is 
significantly  different  from  that  in  Einstein  gravity.  With  this  difference  taken  into  account, 
a  semiclassical  analysis  shows  that  the  oscillating  LQC  solutions  are  still  unstable  with 
respect  to  quantum  collapse.  Our  conclusions  are  summarized  and  discussed  in  Section  V. 


II.  CLASSICAL  SOLUTIONS 

We  consider  a  spatially  hat  but  closed  universe,  which  has  a  toroidal  topology, 


ds2  =  N2(t)dt2  —  a2(t)dx2, 


(1) 
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where  a(t)  is  the  scale  factor  and  N(t)  is  the  lapse  function.  We  shall  use  the  proper  time 
gauge  N  —  1  throughout  the  paper.  Here,  the  Cartesian  coordinates  xl  vary  in  the  range 

0  <  x*  <  L,  (2) 


where  L  =  const  and  the  points  xl  =  0  and  x1  =  L  are  identified.  The  dynamics  of  this 
model  in  LQC  is  accurately  described  by  the  effective  Hamiltonian  [9jj] 

3  cos2  (05). 


n  = 

where  V  =  a3L 3  is  the  volume, 


£2 


87tG72 

0  =  4:irG'ypv, 


V  +  H 


matter  ? 


(3) 


(4) 


and  pv  is  the  canonical  momentum  conjugate  to  V.  Here,  7  is  the  Barbero-Immirzi  parame¬ 
ter,  and  £2  =  A'Ky/Z'yG.  It  is  usually  assumed  that  7  ~  1;  then  £  is  comparable  to  the  Planck 
length. 

We  consider  a  model  where  the  energy  density  is  due  to  a  massless  scalar  held  0  and  a 

negative  cosmological  constant,  A  <  0.  The  corresponding  matter  Hamiltonian  is 
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Umatter  =  ^  +  =  PV,  (5) 

where  is  the  momentum  conjugate  to  0,  and  we  have  defined  the  matter  energy  density 
p.  The  classical  equations  of  motion  for  the  canonical  variables  are  then 


Pp 

<i> 

V 


0 

p± 

V 


3  cos  (0/3) 

7  £ 


sin  (£j8) 


where  in  the  final  relation  we  have  used  the  Hamiltonian  constraint 


(6) 

(7) 

(8) 
(9) 


n  =  o. 


(10) 


Squaring  the  equation  for  V  and  combining  with  Eqs.  ([BJ  and  (10),  we  find  the  modified 
Friedmann  equation: 

87 tG  ( _  p 


a 2 


V 

9H2 


-P 


Pr 


(11) 


1  The  effective  Hamiltonian  commonly  used  in  the  LQC  literature  is  given  by  ([3])  with  cos (£/3)  replaced  by 
sin(03).  These  two  forms  of  H  are  related  by  a  simple  change  of  variable  f)  — >  (3  +  7t/2 l. 
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where 


a/3 


Pmax  327T273 G2  (12) 

is  the  maximum  energy  density  that  can  be  attained  in  this  modclj^]  We  shall  assume  that 

I A I  r$  Pmax- 

The  equations  of  motion  for  0  and  p $  in  Eq.  (|6]{7])  result  in  the  familiar  wave  equation 
for  a  homogeneous  massless  scalar  field, 


a  ■ 

+  3-0  =  0 
a 


Using  Eq.  ([5]) ,  the  matter  energy  density  can  be  expressed  as 


A  C 

P  ~  A+  y2’ 


(13) 


(14) 


where  C  =  p|/2  =  const. 

With  this  matter  density,  the  Friedmann  equation  @  has  an  oscillating  solution  0: 

V(t)  =  Vmin(2X)~1/2  (-  cos(o ;t)  +  1  +  2A)1/2  ,  (15) 

oj  =  \Z96nGpmax\  (1  +  A),  (16) 

where  A  =  |A| / pmax  and 

C 


V min. 


1/2 


=  TSn3 
—  ^  u min’ 


(17) 


v  Pmax  +  |A|_ 

The  universe  oscillates  at  frequency  u  between  minimum  volume  Vmin,  where  p  =  pmax ,  and 
maximum  Vmax  =  Vmin^,  where  p  =  0. 

The  value  of  the  parameter  L  depends  on  the  normalization  of  the  scale  factor  a.  We 
shall  choose  it  so  that  amin  =  1.  Then  Vrnin  =  L3  and  L  =  (C / pmax)1^ {1  +  A)-1/6. 

We  may  also  find  a  solution  for  the  momentum  (3: 


=  (^nGl*Pmax]1/2  arctan 


1  +  A  (  ut\ 

—tan(yj 


We  note  that  in  the  limit  A  — »  0  the  solution  ( 15 )  goes  into 

r(i)  =  dS!>(l  +  24IrGpra««2)I/2, 


(18) 


(19) 


where  V^n  =  {C / Pmax)1^2  ■  This  solution,  which  describes  a  contracting,  bouncing  and 
re-expanding  universe,  has  been  discussed  earlier  by  a  number  of  authors 


2  The  same  form  of  the  modified  Friedmann  equation  is  obtained  in  the  Shtanov-Sahni  braneworld  model 
with  a  timelike  extra  dimension  mum. 
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III.  PERTURBATIVE  STABILITY 


We  generally  expect  to  find  particle  production  in  an  oscillating  universe.  In  our  model, 
this  would  manifest  itself  in  unstable  (growing)  inhomogeneous  fluctuations  of  the  scalar 
field  or  the  metric.  We  shall  now  check  for  this  instability. 

The  dynamics  of  perturbations  in  LQC  has  been  studied  by  Agullo  et  al  pa.  with  the 
conclusion  that  it  is  accurately  given  by  the  usual  quantum  field  theory  on  the  FRW  back¬ 
ground  described  by  the  solutions  to  the  effective  field  equations]^]  In  our  case,  the  scalar 
field  perturbations  should  then  satisfy  the  Klein-Gordon  equation 

a  ■  1 


□0  =  0  +  3-0 - -V  0  =  0 

a  az 


(20) 


in  the  oscillating  background  (15).  Metric  perturbations,  describing  gravitons,  satisfy  the 


same  equation,  so  it  will  be  sufficient  to  investigate  the  stability  of  the  solutions  of  Eq.  (20). 

It  will  be  convenient  to  introduce  a  new  variable  r  =  cot.  The  field  0  can  then  be  expanded 
into  plane  waves, 


0(xw)  =  ^Mr) 

k 

where  the  mode  functions  0*,  satisfy 

k2 


vikx 


+  3 — 00  + 


h  '  9  o 

a  ujzaz 


4>k  =  0. 


(21) 


(22) 


Here,  primes  stand  for  derivatives  with  respect  to  r. 

Periodic  boundary  conditions  in  the  range  (J2])  require  that  /q  =  27rn;/L,  where  nl  are 
integers.  The  eigenvalues  of  the  Laplacian  are  then  given  by 


k  =  (2tt/L)  V, 


(23) 


where  n2  =  n2  +  n2  +  ruj . 


The  scale  factor  corresponding  to  the  solution  ( 15 )  can  be  written  as 

a(r)  =  (2A)-1/6  (—  cos r  +  1  +  2A )1^  , 


(24) 


where  we  have  normalized  to  amin  =  1.  With  this  form  of  the  scale  factor,  the  mode  equation 


(22)  becomes 


+ 


sinr 

27M 


+ 


/1/3(r 


"0fc  0, 


(25) 


3  It  is  shown  in  Ref.  |H]  that  this  prescription  can  become  inaccurate  in  the  regime  where  C  <C  102!6pn 
so  that  the  bounce  is  at  a  near-Planckian  volume  Vmin  <C  10^3.  Here  we  shall  assume  that  C  10 2£6pr] 
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where 


and 


g2  = 


/(t)  =  —  cos  r  +  1  +  2A 
(2A  )V3fc2 


7T 


n 


or 


12  GL2pma;r(2A)2/3(l  +  A)  ’ 


(26) 


(27) 


Eq.  (25)  is  a  form  of  Hill’s  equation,  which  is  a  more  general  form  of  the  familiar  Mathieu 
equation.  As  with  the  Mathieu  equation,  we  can  find  stable  and  unstable  regions  in  the 
parameter  space.  The  modes  </>&  represent  excitations  of  the  field  at  a  certain  momentum  k, 
so  growing  solutions  to  the  mode  equation  signal  particle  production  in  that  state. 

With  a  transformation  to  conformal  time  dr]  =  a~1dt  and  defining  y{rf)  =  a{rf)<p{i]), 


Eq.  (25)  can  be  brought  to  the  form 

y"  +  (k2  +  a"  /  a)y  =  0, 


(28) 


where  primes  stand  for  derivatives  with  respect  to  rj.  One  can  then  use  the  standard  stability 
analysis  using  Floquet  theory  [16].  We  did  not  follow  this  path  because  we  could  not  obtain 
an  analytic  solution  for  a{j])  and  using  a  numerical  solution  would  take  an  excessive  amount 
of  computer  time.  We  therefore  modified  the  Floquet  method  as  described  in  Appendix 


A,  so  that  it  can  be  directly  applied  to  Eq.  (25).  The  resulting  stability  diagram  in  the 
parameter  space  of  q  and  A  is  shown  in  Fig.  1. 

The  parameter  A  in  our  model  determines  the  ratio  amax / amin.  For  A  ~  1,  we  have 
Umax  ~  CLmin,  so  the  scale  factor  oscillates  about  the  value  a  ~  1  with  an  amplitude  also 
~  1.  In  this  regime,  the  diagram  exhibits  the  characteristic  pattern  of  narrow  parametric 
resonance,  with  instability  confined  to  narrow  bands.  For  A  <g7  1,  amax/amjn  ~  A_1//3  ^  1,  so 
the  size  of  the  universe  changes  by  many  orders  of  magnitude  in  the  course  of  one  oscillation. 
The  instability  bands  become  broader  in  this  limit. 

These  results  can  be  qualitatively  understood  as  follows.  In  the  absence  of  oscillations, 
the  energy  spectrum  of  a  scalar  field  in  a  compact  universe  is  discrete.  If  the  oscillation 
amplitude  is  relatively  small,  the  oscillations  act  as  a  periodic  perturbation  of  frequency  c j, 
and  particle  production  occurs  only  if  oj  is  very  close  to  one  of  the  resonant  frequencies  of 
the  modes.  As  the  oscillation  amplitude  gets  large,  the  perturbation  effectively  includes  a 
wide  range  of  frequencies,  so  a  larger  number  of  modes  are  affected. 

As  unstable  mode  functions  oscillate,  their  amplitudes  grow  exponentially, 


4>k{r)  oc  ec 


(29) 
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FIG.  1.  Stability  diagram  in  the  parameter  space  of  q  and  A.  Red  bands  are  unstable. 


with  the  rate  of  growth  a  getting  higher  as  A  gets  smaller  and  the  instability  band  widens. 
The  time  evolution  of  the  energy  density, 

'2  /  k2 


Pk{r)  =  y 


h2  + 


u2a2 


4>k 


(30) 


is  shown  in  Fig.  2  for  q  —  0.55  and  A  =  0.5,  0.05,  0.01.  The  corresponding  growth  rates  are 
a  =  0.09,  0.23,  0.31,  respectively. 


The  parameter  q  defined  in  Eq.  (27)  depends  on  A,  as  well  as  on  L  and  the  mode  number 
n2.  In  Fig.  3  we  show  regions  of  instability  for  the  independent  parameters  of  the  model,  A 
and 

«  =  (.GL2pmax)~1/2.  (31) 


Each  unstable  band  of  Fig.  1  now  splits  into  an  infinite  number  of  bands,  corresponding  to 
different  values  of  n2.  We  have  included  in  Fig.  3  only  modes  with  n2  <  10,  with  higher 
values  of  n2  indicated  by  lighter  shades  of  grey. 

It  is  interesting  to  compare  our  results  with  those  of  Graham  et  al  [5],  who  studied  the 
stability  of  their  simple  harmonic  universe  model.  In  that  model,  the  solution  for  a(r)  is 
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A=.01,. 05,.5,  q=.5 

Inpft(r) 


FIG.  2.  The  growth  of  energy  density  in  unstable  modes  with  q  =  0.55  and  A  =  0.5,  0.05,  0.01. 
Lower  curves  on  the  graph  correspond  to  larger  values  of  A. 


inversely  proportional  to  the  oscillation  frequency  oj,  so  u>  drops  out  of  the  mode  equation 


(25).  This  leaves  a  single  model  parameter  7  =  3 1  A| / 27rG'p§,  where  A  <  0  is  the  cosmological 
constant  and  po  characterizes  the  contribution  of  matter  with  w  =  —2/3  to  the  total  density. 
The  role  of  this  parameter  is  similar  to  that  of  A  in  our  model.  Graham  et  al  find  that  for 
7  ~  1,  amax/amin  ~  1  and  all  modes  are  stable,  while  for  7  <C  1,  amax/amin  1  and 
there  is  a  large  number  of  unstable  modes.  (An  oscillating  solution  exists  only  for  7  <  1.) 
In  contrast,  our  diagram  in  Fig.  3  shows  a  non-trivial  pattern  of  stability  and  instability 
regions.  In  particular,  we  find  wide  ranges  of  k  where  the  model  is  unstable  when  A  ~  1 
and  where  it  is  stable  when  A  is  very  small. 


IV.  TUNNELING  TO  NOTHING 

Even  if  the  oscillating  LQC  universe  is  stable  with  respect  to  classical  perturbations,  it 
may  tunnel  quantum  mechanically  from  the  bounce  point  at  V  =  Knm  to  a  state  of  zero 
volume.  Here  we  shall  estimate  the  semiclassical  tunneling  probability. 

The  application  of  standard  semiclassical  methods  to  this  case  is  complicated  by  the  fact 
that  the  held  0  has  a  nonzero  velocity  at  the  bounce  point,  0(14im)  =  Pp/Knm  [see  Eqs.  (Jg]- 
0]-  This  problem,  however,  is  not  difficult  to  resolve  if  we  note  that  the  variables  V  and  0 
are  essentially  decoupled  from  one  another.  The  Hamiltonian  ([3]),(|5|  is  independent  of  0,  so 
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FIG.  3.  Stability  diagram  in  the  parameter  space  of  k  and  A,  with  instability  regions  corresponding 
to  different  values  of  n2  marked  by  different  shades  of  grey.  Lighter  shades  correspond  to  larger 
values  of  n2. 


Ptf,  =  (2 C)1//2  =  const,  and  the  wave  function  can  be  factorized 

T(V,0)  =  eip^{V).  (32) 

The  problem  then  reduces  to  a  one- dimensional  tunneling  problem  with  a  single  dynamical 
variable  V,  which  is  described  by  the  action 

s  =  I'dt{^v-Nn)’  <33> 

where 

n  =  -Pmaxv  cos2  m  +  £  +  AV.  (34) 

The  semiclassical  tunneling  probability  is  given  by 

V  ~  e~2SE,  (35) 


where  Se  is  the  Euclidean  action  of  the  instanton  solution  to  the  Euclidean  equations  of 


motion.  The  instanton  can  be  obtained  from  the  Lorentzian  solution  (15), (18)  by  analytic 
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continuation  t  — >  —it,  (3  — >  i/3  to  the  classically  forbidden  range  of  0  <  V  <  Vmin.  This  gives 


V{t)  =  L3( 2A)"1/2  (-  cosh (ut)  +  1  +  2A)1/2  , 


jd(t)  =  -  (87rG7*PmaxY/2  tanh~ 1 


1  +  A  /  cot 

— tanh  y 


(36) 

(37) 


The  classically  forbidden  range  extends  from  t  =  0,  where  f3  =  0,  to  t  =  tf,  where  V  =  0. 
( tf  can  be  found  from  cosh(cnf/)  =  1  +  2 A.)  Taking  into  account  the  Hamiltonian  constraint 
T-L  =  0,  the  instanton  action  is  given  by 


Sk  = 


rtf 


AttG^ 


dt 


pdZ 

dt 


A'lrG'y 



dV/3(V) 


(38) 


The  form  of  /3{V)  may  be  determined  from  the  constraint, 

${V)  =  r 1  cosh'1 


C 


Pm£\!  ~  lAl 


(39) 


Then,  introducing  a  new  variable  x  =  V/\ 3min,  the  integral  (38)  can  be  rewritten  as 

1/2 


SF  = 


14 


AitG^/i 


dx  cosh  1  ( — —  A 
x 2 


i  o 


(40) 


It  can  be  expressed  as  an  elliptic  integral,  but  this  expression  is  not  particularly  illuminating 
and  we  do  not  present  it  here.  An  interesting  special  case  is  the  limit  A  — »  0,  when  Vmax 


Vmin.  In  this  limit  we  obtain 


V  ~  exp 


71 


a/3  14 


2  P 


where  we  have  used 


dx  cosh  1  —  =  — . 

x  2 


(41) 


(42) 


The  above  semiclassical  treatment  indicates  that  when  the  universe  reaches  its  minimum 
volume,  there  is  a  non-zero  probability  of  tunneling  to  a  singularity.  The  tunneling  is 
strongly  suppressed  when  the  minimum  volume  at  the  bounce  14im  is  much  larger  than  the 
Planck  volume. 

We  note  finally  that  our  conclusions  here  are  somewhat  different  from  those  of  Ashtekar  et 
aim  ,  who  studied  quantum  tunneling  in  the  same  LQC  model  with  A  =  0.  They  found  that 
the  Euclidean  action  in  the  classically  forbidden  region  between  V  =  0  and  V  =  Vmin  is  Se  — 
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0,  suggesting  that  the  tunneling  to  V  =  0  is  unsuppressed.  On  the  other  hand,  numerical 
calculations  in  Ref.  B33  indicate  that  the  wave  function  is  actually  suppressed  in  this  region]/] 
and  the  authors  interpret  this  as  a  breakdown  of  the  semiclassical  approximation.  In  our 
view,  the  semiclassical  approximation  is  accurate  under  the  usual  conditions  (roughly,  Se 
1).  The  reason  for  the  discrepancy  is  that  the  Euclidean  continuation  in  [T7]  was  performed 
in  the  full  action,  including  both  V  and  (ft  variables,  while  we  considered  a  reduced  action 


(|33|)  with  a  single  variable  V.  The  latter  appears  to  be  the  correct  prescription  in  the 
presence  of  classical  motion  in  A'^ 


V.  SUMMARY  AND  DISCUSSION 

A.  Particle  production 

We  investigated  the  stability  of  some  oscillating  universe  solutions  to  the  effective  equa¬ 
tions  in  Loop  Quantum  Cosmology.  These  solutions  have  the  topology  of  a  3-torus  and 
contain  matter  in  the  form  of  a  homogeneous  massless  scalar  held  and  a  cosmological  con¬ 
stant  A  <  0.  We  found  that  for  some  values  of  the  parameters  the  model  is  stable  with 
respect  to  small  perturbations  of  the  scalar  held,  indicating  that  oscillation  of  the  scale  fac¬ 
tor  does  not  result  in  particle  production.  There  are,  however,  substantial  regions  of  the 
parameter  space  which  are  unstable. 

In  a  universe  described  by  an  unstable  solution,  there  will  be  particle  production.  The 
energy  density  of  the  created  particles  will  grow  exponentially  with  time  and  will  eventually 
become  comparable  to  that  of  the  homogeneous  held.  Our  perturbative  analysis  cannot  be 
extended  beyond  this  point;  here  we  shall  attempt  to  give  a  qualitative  outline  of  subsequent 
evolution.  For  definiteness  we  shall  assume  that  A  1,  so  that  Vmax/Vmin  1. 

The  additional  matter  component  due  to  the  created  particles  can  be  approximated  as  a 
perfect  fluid  with  a  radiation  equation  of  state  w  =  1/3,  so  the  total  energy  density  is 

P  =  A  +  —  +  ^^3  =  A  +  p<p  +  pr,  (43) 

where  D(t)  grows  with  time.  As  the  universe  contracts,  the  ‘radiation’  energy  density  pr 
grows  slower  than  p^,  so  it  will  have  little  effect  on  the  dynamics  of  the  bounce  at  «mm,  at 

4  We  also  mention  a  related  result  by  Craig  [T5J  who  showed  analytically  that  eigenfunctions  of  the  quantum 

evolution  operator  in  LQC  decay  exponentially  in  the  region  between  zero  volume  and  the  bounce. 

5  For  a  discussion  of  multi-dimensional  tunneling  in  the  presence  of  classical  motion,  see  Ref.  [19]. 
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least  initially.  By  the  same  token,  pr  decreases  slower  than  p ^  during  the  expansion  and  will 
eventually  exceed  p $  at  large  a.  As  a  result,  it  will  take  longer  to  reach  the  bounce  at  amax 
(where  p  =  0),  and  the  value  of  amax  will  increase.  The  effect  of  this  is  similar  to  that  of 
decreasing  A. 

If  the  initial  state  at  the  onset  of  instability  is  in  one  of  the  narrow  bands  in  Fig.  3,  then 
any  change  of  the  parameters  is  likely  to  take  us  out  of  the  resonance  and  to  stabilize  the 
oscillating  universe.  On  the  other  hand,  if  the  initial  state  is  in  a  broad  instability  band  (say, 
at  a  relatively  small  k  ~  0.05  —  0.1),  then  an  effective  decrease  of  A  can  make  the  solution 
even  more  unstable.  Eventually  pr  will  become  larger  than  p $  even  at  small  a  ~  aTOjn.  The 
bounce  will  then  occur  at  a  larger  value  of  Vmm ,  which  amounts  to  an  effective  increase  in 
L  and  a  decrease  in  k.  A  look  at  the  diagram  in  Fig.  3  suggests  that  if  we  start  at  relatively 
small  values  of  k  and  A  and  move  in  the  direction  of  still  smaller  n  and  A,  we  may  never 
leave  the  instability  region.  This  means  that  particle  creation  may  be  a  runaway  process, 
with  both  amax  and  amin  becoming  larger  and  larger  and  more  and  more  particles  being 
produced.  Whether  or  not  this  happens  depends  on  the  initial  state,  the  other  option  being 
a  stable  oscillating  universe.  This,  however,  is  not  the  whole  story. 


B.  Quantum  collapse 

Even  if  the  oscillating  universe  is  stable  with  respect  to  particle  production,  we  found 
that  every  time  it  bounces  at  a  =  amin,  there  is  a  nonzero  probability  for  it  to  tunnel  to  zero 
size,  a  =  0.  Classical  solutions  of  GR  become  singular  at  a  =  0,  and  it  is  usually  assumed 
that  the  spacetime  terminates  there.  Then  a  tunneling  to  a  =  0  would  be  the  end  of  the 
universe. 

An  alternative  scenario  is  that  the  singularity  at  a  =  0  will  eventually  be  resolved  in  the 
quantum  theory  of  gravity.  It  could  be  replaced,  for  example,  by  a  Planck-size  nonsingular 
nugget  with  a  ~  t.  Then  a  tunneling  am*n  — >  l  could  be  followed  by  the  inverse  process 
t  — >  amin,  and  the  oscillation  would  resume.  If  this  state  of  affairs  can  continue  indefinitely, 
this  could  serve  as  a  basis  for  a  viable  model  of  an  emergent  universe. 

However,  if  a  tunneling  t  — >  amin  is  possible,  it  seems  likely  that  the  parameters  C 
and  D  characterizing  the  energy  content  of  the  resulting  oscillating  universe  will  generally 
be  different  from  their  initial  values  at  the  time  of  quantum  collapse.  The  cosmological 
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constant  A  may  also  be  different  if  the  particle  physics  model  admits  a  landscape  of  vacua 
with  different  values  of  A.  The  universe  will  then  explore  the  entire  parameter  space  and 
will  eventually  get  into  the  regime  of  runaway  particle  production. 

The  following  argument  suggests  that  quantum  production  of  increasingly  large  numbers 
of  particles  may  be  a  general  feature  of  any  model  of  a  compact  nonsingular  universe.  The 
conserved  energy-momentum  ‘complex’  Qv  in  GR  can  be  expressed  in  terms  of  a  ‘superpo- 
tential’  U"a  = -U™  (20], 

%  =  daUf,  (44) 

dv%  =  0.  (45) 

The  total  energy  and  momentum  of  a  compact  universe  are  then  identically  zero, 

P,,  =  J  d3xdiU f  =  0,  (46) 

with  negative  gravitational  energy  exactly  compensating  the  energy  of  matter.  Hence,  the 
production  or  destruction  of  particle-antiparticle  pairs  in  such  a  universe  are  not  forbidden 
by  any  conservation  law.  As  a  rule  of  thumb  in  quantum  mechanics,  any  process  that  is  not 
strictly  forbidden  will  occur  with  some  probability.  Furthermore,  an  increase  in  the  number 
of  particles  increases  the  entropy  of  the  universe.  This  indicates  that  a  compact  universe 
will  evolve  towards  states  with  a  larger  number  of  particles.  In  oscillating  models  of  the  kind 
we  considered  here,  this  is  accompanied  by  an  increase  in  the  size  of  the  universe.  Given  an 
infinite  time,  the  number  of  particles  and  the  size  of  the  universe  would  then  grow  without 
bound.  A  quantitative  analysis  of  these  issues  should  await  progress  in  the  quantum  theory 
of  gravity. 
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APPENDIX 


In  this  Appendix  we  shall  derive  a  stability  criterion  for  solutions  to  the  equation 


a  ■  k2 

4>k  +  3— H — ~4>k  =  0, 
a  a - 


(47) 


where  the  scale  factor  a(f )  is  a  periodic  function  with  period  T. 


According  to  Floquet’s  theorem,  Eq.  (47)  admits  solutions  of  the  form 


<p(t)  =  e“p(«), 


(48) 


where  p{t)  is  a  periodic  function  with  period  T,  and  a  is  a  constant  defined  by  the  boundary 
conditions.  It  is  apparent  that  when  a  is  real,  the  solutions  are  oscillatory  and  therefore 
stable,  whereas  when  a  is  complex,  the  solutions  grow  or  decay  and  are  unstable. 


Note  that  if  0(f)  is  a  solution  of  (47),  then  0(f  +  T)  is  also  a  solution.  With  the  ansatz 
I),  we  have 

0(f  +  T)  =  C0(f),  (49) 

where  (  =  eiaT  is  a  constant. 

To  derive  the  stability  criterion,  we  first  define  two  solutions,  0i(f)  and  02(f),  by  the 
initial  conditions 

0r(O)  =  1  02(O)  =  0 

0i(O)  =  0  02(0)  =  1.  (50) 

Any  solution  can  be  written  as  a  linear  combination  of  these  two  solutions: 


0(f)  =  Cl  01  (f)  +  C202  (f)  • 

In  particular,  the  solutions  0i(f  +  T)  and  02(f  +  T)  are 

0i(f  +  T)  =  0i(T)0i(t)  +  0i(T)02(f) 

<h(t  +  T)  =  <p2(T)<pi(t)  +  <p'2(T)<p2(t). 


(51) 


(52) 


For  any  solution  of  the  form  (49),  we  obtain,  using  Eqs.  (51),  (52),  the  following  set  of 
linear  equations  for  Ci  and  c2: 


(0i(T)-C)ci  +  02(T)c2  =  O 
01(T)ci  +  (0'(T)-C)c2  =  O 


(53) 

(54) 
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Nonzero  solutions  to  this  set  of  equations  exist  when 


=  C2  -  (MT)  +  0'2(T))C  +  W(MT),  02 (T ))  =  0  (55) 


(MT)-C)  MT) 

MT)  (MT)-  0 

where  W(t)  =  W(0i,  (t),  02(£))  is  the  Wronskian.  From  Eq.  ([47]),  W(t)a3(t)  =  const,  so 


that  W(T)  =  W( 0).  With  the  boundary  conditions  in  Eq.  (50),  we  have  W( 0)  =  1,  so  that 
Eq.  (|55j)  becomes 

c2  -  (MT)  +  MT))(  +  1  =  0.  (56) 

The  product  of  the  two  roots  of  this  equation  is  equal  to  1;  hence  the  roots  can  be 
represented  as  £1,2  =  exp (±iaT).  The  sum  of  the  roots  is 

Cl  +  C2  =  2  cos(aT)  =  MT)  +  MT)  =  b.  (57) 


If  b  >  2,  then  the  roots  are  real,  while  if  b  <  2,  the  roots  are  complex  and  conjugate  to  one 
another. 

As  we  indicated  above,  the  solutions  are  stable  for  real  a  and  unstable  for  complex  a. 
This  corresponds  to  the  condition 

|01(T)  +  0,2(T)|<2  (58) 


for  stable  solutions,  and 

|01(T)  +  C4(T)|>2  (59) 

for  unstable  solutions^  In  order  to  check  the  mode  stability  in  our  model,  we  integrate 


Eq.  (47)  with  a(t)  from  Eq.  (24),  using  the  initial  conditions  in  (50).  Then  we  check  whether 


the  solutions  at  t  =  T  satisfy  (58)  or  (59). 

The  stability  diagram  in  Fig.  3  was  produced  by  sampling  the  parameter  space  0  <  A,  n  < 
1  at  logarithmic  intervals,  so  there  are  many  more  points  at  small  values  of  the  parameters 
than  at  A,  n  ~  1.  This  ensures  that  we  have  high  resolution  in  regions  where  it  is  required. 
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